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Report  Summary 


One  of  the  most  important  physical  phenomena  is  resonance.  Musical  instru¬ 
ments,  radio,  and  television  depend  on  tuning  strings,  pipes,  or  circuits  to  resonance 
at  a  particular  frequency.  At  resonance  the  amplitude  of  vibrations  or  oscillations 
reaches  a  maximum  due  to  the  progressive  amplification  of  multiple,  in-ph^ise  reflec¬ 
tions.  The  purpose  of  Contract  F19628-88-K-0024  has  been  to  explore  the  properties 
and  possible  applications  of  a  remarkable  new  resonator  that  combines  the  charac¬ 
teristics  of  an  antenna  array  with  those  of  a  cavity  resonator.  It  consists  of  a  large 
number  of  parallel  coplanar  dipoles,  equally  spaced  along  a  closed  loop  that  may  be 
circular,  elliptical  or  serpentine.  The  single  driven  element  excites  waves  of  electric 
current  in  the  elements  that  travel  in  both  directions  around  the  loop.  The  loop  can 
be  tuned  to  resonance  by  correctly  choosing  the  length  and  radius  of  each  element, 
the  distance  between  them,  and  their  total  number.  When  the  loop  is  resonant, 
extremely  large  currents  oscillate  in  all  of  the  elements  with  phase  relations  char¬ 
acteristic  of  any  one  of  a  large  number  of  possible  modes.  These  generate  a  field 
pattern  determined  by  the  phase  relations  of  the  currents  and  the  particular  shape  of 
the  closed  loop.  Unlike  wire  and  waveguide  resonators  in  which  the  currents  travel 
the  entire  length  of  the  circuit  with  correspondingly  large  losses,  the  open  resonator 
involves  only  transverse  currents  in  very  short,  widely  spaced  conductors.  Possible 
applications  as  a  resonant  circuit  include  oscillators  and  frequency  control  where 
its  unusual  shape  and  extremely  low  losses  make  it  useful.  Since  the  conditions  of 
resonance  are  largely  independent  of  the  shape  of  the  closed  loop  so  long  as  there  is 
no  radius  of  curvature  smaller  than  two  wavelengths,  an  enormous  range  of  possi¬ 
ble  shapes  and  associated  field  patterns  is  available.  These  may  include  particular 
pheise  relations,  shapes  and  sizes  that  lead  to  extremely  high  directivity.  Such  an 
array  would  have  !"^poitant  practical  applications  in  point-to-point  communication 
and  ru,Jar.  Another  application  is  as  a  niicrowave  beacon. 

The  technical  problems  in  the  research  are  theoretical  and  experimental.  They 

iii 


are  related  to  extremely  precise  design  and  construction.  The  detailed  mathemati¬ 
cal  analysis  combined  with  quadruple-precision  numerical  calculation  has  provided 
designs  of  which  at  least  one  has  been  verified  experimentally  (supported  by  another 
contract).  The  general  methodology  involves  a  highly  precise  analytical  formula¬ 
tion,  its  accurate  evaluation  by  computer  and,  finally,  a  meticulous  verification  using 
an  experimental  model  constructed  to  high  dimensional  tolerances.  The  results  ob¬ 
tained  to  date  have  been  highly  encouraging.  They  show  conclusively  that  this 
novel  open  resonator  can  be  designed  and  built  and  that  it  is  a  device  with  very 
unusual  properties  that  can  have  applications  of  major  importance. 

The  design  of  resonant  arrays  in  circular  form  has  been  achieved.  The  problem 
of  further  refinement  to  sharpen  the  resonances  remains  to  be  studied,  including 
the  possible  need  for  superconducting  elements.  Of  primary  interest  for  future 
research  is  the  systematic  study  of  the  effects  of  changes  in  shape  on  the  field 
pattern.  Can  a  superdirective  array  be  achieved?  The  unique  properties  of  single¬ 
element  excitation,  very  large  resonant  currents,  purely  resistive  input  impedance, 
and  extremely  sharp  pencil-like  beams  already  observed  for  the  circular  array  are 
encouraging. 
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I.  Review  of  Work  Completed  During  the  Period 
March  28,  1988  to  September  30,  1990 

The  study  of  the  properties  of  closed  loops  of  parallel  coplanar  dipoles  be¬ 
gan  with  a  systematic  review  of  early  work  on  superdirective  endfire  arrays,  long 
parasitic  arrays  of  Yagi  type,  and  the  resonemt  properties  of  circular  arrays  with 
only  one  element  driven.  This  summarizing  investigation  is  contained  in  the  paper 
“Supergain  Antennas  and  the  Yagi  and  Circular  Arrays”  by  R.  W.  P.  King  [1].  It 
revealed  that  there  is  a  minimum  attenuation  design  for  Yagi-type  structures  with 
a  hundred  or  more  elements  and  that  this  applies  not  only  to  linear  arrays  but  to 
structures  of  semicircular  form  provided  the  radius  of  curvature  is  large  enough. 
When  two  semicircular  arrays  are  combined  to  form  a  circular  array  driven  by  a 
single  element,  traveling  waves  of  current  can  be  excited  that  propagate  around  the 
circle  in  opposite  directions.  Clearly,  with  proper  dimensions,  the  conditions  for 
a  resonant  system  are  provided  in  which  the  currents  successively  induced  in  the 
elements  can  build  up  very  large  amplitudes. 

It  W£LS  initially  thought  that  elements  near  an  antiresonant  length  were  appro¬ 
priate  and  a  systematic  study  of  the  three-term  theory  was  undertaken  as  a  tool 
for  analyzing  circular  arrays  with  elements  up  to  a  wavelength  long.  In  the  first 
part  of  this  investigation  the  properties  of  the  vector  potentials  of  coupled  anten¬ 
nas  were  examined  and  used  to  show  that  great  simplification  is  achieved  when  the 
distance  between  adjacent  elements  is  equal  to  or  greater  than  the  half-length  of 
the  elements.  This  is  contained  in  the  paper  “Electric  Fields  and  Vector  Potentials 
of  Thin  Cylindrical  Antennas”  by  R.  W.  P.  King  [2].  The  complete  analysis  of 
circular  arrays  in  which  the  three- term  theory  is  applied  to  determine  the  currents 
in  the  elements  is  contained  in  the  paper  “The  Large  Circular  Array  with  One  Ele¬ 
ment  Driven”  by  R.  W.  P.  King  [3].  An  important  conclusion  was  that  the  simpler 
two-term  theory  is  entirely  adequate  when  the  elements  are  electrically  short. 
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Since  the  trigonometric  functions  used  in  the  two-  and  three-term  theories  can¬ 
not  take  account  of  the  end  effect  that  characterizes  actual  currents  on  tubular 
elements,  a  detailed  analysis  was  made  to  improve  the  currents  near  the  open  ends. 
This  included  the  derivation  of  “The  Universal  Current  Distribution  near  the  End 
of  a  Tubular  Antenna”  in  a  paper  by  H.-M.  Shen  and  T.  T.  Wu  [4]  and  “The  Com¬ 
bination  of  the  Universal  End-Current  and  the  Three- Term  Current  on  a  Tubular 
Antenna”  in  a  paper  by  H.-M.  Shen,  R.  W.  P.  King,  and  T.  T.  Wu  [5]. 

The  application  of  the  two-term  theory  to  60-  and  90-element  circular  arrays 
of  electrically  long  and  thin  elements  yielded  no  sharply  resonant  arrays.  On  the 
other  hand,  the  two-term  theory  applied  to  60-  and  90-element  circular  arrays  of 
electrically  short  and  quite  thick  tubular  elements  led  to  the  discovery  of  a  variety 
of  sharply  resonant  modes  with  large  currents  in  all  elements  and  radiation  patterns 
consisting  of  sharp  spikes.  These  modes  and  the  underlying  theory  are  described  in 
detail  in  the  paper  “The  Resonant  Circular  Array  of  Electrically  Short  Elements” 
by  G.  Fikioris,  R.  W.  P.  King,  and  T.  T.  Wu  [6].  The  properties  of  the  different 
resonant  modes  are  related  to  the  resonances  of  particular  phase  sequences  in  the 
method  of  symmetrical  components  used  in  the  analysis.  They  depend  critically  on 
the  radius  and  length  of  each  element  and  the  distance  between  adjacent  elements. 

A  comparison  of  the  predicted  resonance  in  and  the  associated  field  pattern 
of  a  90-element  circular  array  with  those  observed  in  an  experimental  array  of 
90  brass  monopoles  on  an  aluminum  ground  screen  and  driven  by  one  element 
showed  quite  good  agreement.  However,  the  extremely  sharp  resonances  needed  for 
a  possible  superdirective  array  clearly  require  great  refinement  in  both  the  theory 
and  experiment  in  order  to  achieve  much  higher  accuracy.  Recent  very  significant 
progress  in  the  theoretical  treatment  is  described  in  the  following  section. 
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II.  Recent  Progress  in  the  Theoretical  Analysis  of 
Resonant  Loops  of  Parallel,  Coplanar  Tubular  Dipoles 

Some  aspects  of  current  work  are  summarized  in  the  paper  “Analytical  Studies 
of  Large  Closed-Loop  Arrays”  by  G.  Fikioris,  D.  K.  Freeman,  R.  W.  P.  King,  H.- 
M.  Shen,  and  T.  T.  Wu  [7],  which  was  presented  at  the  recent  SPIE  0E/LASE’91 
Meeting  held  in  Los  Angeles,  January  20-25,  1991.  Most  of  this  material  is  included 
in  the  complete  review  of  recent  and  current  work  given  in  Section  III. 

An  important  feature  of  the  improved  analysis  is  that  a  way  has  been  found 
for  relating  an  infinite  mmiber  of  integral  equations  to  a  relatively  small  number. 
This  permits  great  simplification  in  the  analysis  of  the  resonant  circular  array  since 
the  kernel  of  the  finite  number  of  integral  equations  can  be  approximated  by  the 
much  simpler  kernel  of  the  infinite  linear  curay.  Specifically,  the  limit  of  the  kernel 
is  the  kernel  of  the  infinite  linear  array;  indeed  it  should  be  expected  that  the 
solution  for  the  circular  array  with  large  radius  of  curvature  and  many  elements 
should  be  close  to  that  of  the  infinite  linear  array.  Physically,  the  circular  array 
supports  waves  traveling  around  it  continuously,  in  both  directions  and  with  little 
attenuation.  In  the  infinite  linear  array,  the  so-called  surface  waves  [8]  travel  to 
infinity.  This  establishes  a  definite  connection  between  the  infinite  linear  array — 
which  is  governed  by  an  infinite  number  of  integral  equations — and  the  large  circular 
array — which  is  governed  by  a  finite  number  of  integral  equations. 
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III.  Improved  Analysis  of  the  Resonant  Circular  Array 
of  Electrically  Short  Elements 


A.  Introduction 

In  the  recent  paper  by  Fikioys,  King  and  Wu  [6],  it  was  shown  theoretically 
that  a  large  circular  array  of  electrically  short  dipoles  may  possess  very  narrow 
resonances.  The  array  consists  of  N  identical,  perfectly  conducting,  parallel,  non- 
staggered  elements.  Only  element  1  is  center-driven  with  a  voltage  Kj;  the  other 
N  —  1  elements  are  parasitic.  The  half-lengths  and  radii  of  the  dipoles  are  h  and 
a,  respectively;  it  is  assumed  that  fca  ■C  1  and  a  -C  <  A/2,  where  A  =  27r/fc 
is  the  free-space  wavelength.  The  distance  between  adjacent  elements  is  d  >  h. 
The  existence  of  narrow  resonances  has  been  recently  verified  experimentally  [9]. 
Resonances  occur  when  N  and  the  radius  p  =  d/[2  sin(7r/iV)]  of  the  circle  are  large. 

The  N  currents  satisfy  N  coupled  integral  equations  and  the  approximate  so¬ 
lution  used  in  [6]  W2is  the  two-term  theory  developed  by  Mack  [10],  [11].  Thus, 
the  current  on  each  parasitic  element  is  proportional  to  {cos  kz  —  cos  kh)  and  the 
current  on  the  driven  element  has  an  additional  term  sin  k{h  —  jzj),  where  z  is  the 
distance  from  the  center  of  the  dipole.  In  this  section  it  is  assumed,  for  simplicity, 
that  N  is  even.  It  follows  that  the  N  currents  are  a  superposition  of  the  N/2  -|-  1 
phase-sequence  currents  7^ ’"^(2),  m  =  0,  1,  . . . ,  N/2,  where  is  the  current 

on  element  1  when  all  elements  are  dr’ven  by  voltages 

The  array  is  at  its  phase-sequence  resonance  when  becomes  large. 

Recent  investigations  [12]  show  that  the  kernel  of  the  phase-sequence  inte¬ 
gral  equation  used  in  [6]  does  not  take  into  account  the  coupling  effects  with  enough 
accuracy  for  a  resonant  circular  array.  A  modified,  very  simple  kernel  has  been  de¬ 
veloped  [12];  it  is  obtained  by  setting  a  =  0  in  the  self-term  of  the  imaginary  part 
of  the  kernel  used  in  [6].  This  modified  kernel  is,  under  suitable  conditions,  expo¬ 
nentially  small  in  N,  suggesting  the  existence  of  resonances  that  are  exponentially 
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narrow.  Here  use  is  made  of  the  two-term  theory  developed  by  Mack  and  applied  in 
[6]  but  with  the  modified  kernel  to  analyze  the  circular  array  at  (or  near)  resonance. 
Conditions  for  resonance  are  developed,  the  significance  of  the  many  parameters  in¬ 
volved  is  examined,  and  simple  formulas  for  the  radiation  field  are  derived.  Ways 
to  optimize  performance  are  indicated.  The  resonant  circular  array  is  shown  to  be 
a  highly  directive  structure  with  field  patterns  that  involve  many  sharp,  pencil-like 
beams. 

B.  Two- Term  Theory  and  the  Modified  Kernel 

The  two-term  theory  solution  [6],  [11]  for  the  current  on  element  I  is  valid 
subject  to  the  conditions  noted  in  the  previous  secti'ai.  The  formulas  are  repeated 
here,  with  a  notation  more  convenient  for  the  present  purposes: 


=  { 


j27rVy 


sin  k(h  —  1^1)  -f  Tj(cos  —  cos  kh)];  1  =  1, 


r^(cos  kz  —  cos  kh);  /  =  2,  1,  . . . ,  N, 


(1) 


where  (q  =  —  376.73  Q.  The  parameter  real  and  independent  of 

N  and  d/X.  It  is  defined  in  [6],  [llj.  The  coefficients  T,  of  the  shifted-cosine  part 
of  the  current  are  complex  and  depend  on  all  the  parameters  of  the  problem.  They 
are  obtained  by  superimposing  the  phase-sequence  coefficients 


T, 

where 


N/2-1 

j(o)  _  (_i)'r(^/2)  ^  2  ^  r<'"^cos 

rn  =  l 


2n{l  —  l)m] 
N 


}• 


(2) 


'J'(fn)  _ 


D^n'+jO'; 


(m) 


P  -L 


m  =  0,  1,  . . . ,  N/2. 


(3) 


In  (3),  the  parameters  in  the  numerator  and  the  denominator  are  all  real.  The 
subscript  I  means  that  the  parameter  depends  only  on  the  in; aginary  part  of  the 
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modified  phase-sequence  kernel.  Hence,  and  are  independent  of  the 

radius  a/ A.  The  subscript  IP  means  that  the  parameter  depends  only  on  the  self¬ 
term  of  the  real  part  of  the  kernel  and  is  therefore  independent  of  N,  m  and  d/\. 
The  subscript  SP  means  that  the  quantity  depends  only  on  the  mutual  terms  of 
the  real  part  of  the  kernel  and  is  therefore  independent  of  a/X.  The  full  formulas 
for  the  various  parameters  are 


=  J  ^smk(h  -  -  z)dz. 


X  [cos  A’/i  A'S(2)  -  K^ii\h  -  z)]  dz, 


(”j)/ 


X  [cos  AA  Ki”*\z)  —  A'/’"^(h  —  z)]  dz, 

1  /■'■ 
i-cosk 


X  [cos  kh  K^fi{z)  —  Kinih  —  s)]  dz, 


=  - —  f  (cos  kz  —  cos  kh) 

“  1  —  cos  kh  J_i^ 

:<  [cos  kh  -  ld"ii\h  -z)]  dz, 

=  - —  I  (cos  kz  —  cos  kh) 

1— cos  kh  J_f^ 

X  [cosA,A/v:!'"^(^)  -  ld'"\h  - 


(4) 

(5) 

(6) 

(7) 

(8) 

(9) 


The  various  parts  of  the  modified  kernel  are 


cos  kR^{z) 


N/2+l 

1=2 


'2~{l  —  1  )in 

N 


cos  kRi(z) 
Hiiz)  ■ 


(10) 

(11) 
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.  ,  A72+1 


cos 


1=2 


2n{l  —  l)m 
N 


sin  kRi(z) 


where 


(12) 


?,= 


and 


/  =  N/2  +  1, 
otherwise. 


R,(z)=(z^  +  b-‘)'lh 


a. 


b/ =  ^  <isin[(/ —  l)7r/A^] 
sin(7r/A^) 


/=1, 

1^1. 


(13) 


(14) 


Note  that  the  radius  a  does  not  appear  in  (12).  Finally,  the  relation  between 
and  the  phase-sequence  admittances  is 


y(».)  ^  G^'n)  ^  Jll.  ^  [sin  kh  -f  r<"»>(l  -  cos  kh)],  (15) 

and  the  self-  and  mutual  conductances  G^  ,  (susceptances  By  ,)  are  determined  only 
by  the  phase-sequence  conductances  (susceptances  2?^"*^): 


=  lJr(0)_(_i)/y(/v/2)^2  Y.  r<'">cos 

m=l 

The  modified  kernel  (10)- (12)  has  been  evaluated  asymptotically  subject  to 
the  conditions  N  —*  oo  and  d/X  fixed  (to  be  more  precise,  the  arc-length  spacing 
kp/N  is  fixed)  with  d/\  <  mfN  <  |  [12].  K^l^iz)  is  well  approximated  by  the 
kernel  of  the  infinite  linear  aimy  (replace  by  by  Id  in  (11)  and  let  — >  oo  while 
keeping  rti/N  fixed).  It  is  therefore  roughly  independent  of  N  for  large  N  and  fixed 
m/N .  The  imaginary  part  7\.'J"'^(z)  must  be  approximated  more  carefully  because 
the  imaginary  part  of  the  kernel  of  the  infinite  linear  array  is  exactly  zero.  The 
asymptotic  formula  for  K)  {z)  is 

K\'"\z)  1  1  1 

It  ~  47r'/2  A^»/2  [(m/iV)2 -(d/A)2]3/4 


27r(/-  l)m 

N 


} 


(16) 


exp(-2A'(m/iV)j(i„)) 
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(17) 


+  {same  with  m  N  —  m}, 

where 

~  m/^’  0<x<l  and  <7(3-)  =  cosh~' ^  -  (1  -  (IS) 

The  approximation  is  better  when  z/A  is  small  and  when  d/\  is  not  very  close  to 
m/N.  These  are  the  cases  of  interest.  When  m  <C  iV/2,  only  the  first  term  needs  to 
be  kept;  in  the  extreme  case  of  m  =  iV/2,  the  second  term  simply  contributes  a  factor 
of  2.  The  function  g{x)  appearing  in  the  exponential  is  positive  and  decreasing,  with 
5'(0)  =  00.  The  following  properties  of  are  noted: 

Property  (1):  K)  ^(2)  is  approximately  independent  of  2  when 
d/A  <  m/N. 

Property  (2):  k['^\z)  is  exponentially  small  in  N  for  fixed  cl/X 

and  fixed  m/N  with  d/X  <  m/N. 

Property  (3):  K\^\z)  is  a  rapidly  decreasing  function  oi  d/X 

when  N  and  m/N  are  fixed  with  d/X  <  m/N . 

With  the  use  of  Property  (1),  approximations  for  and  are  obtained 

when  h/X  is  not  too  large.  Thus, 

_  _2(sin  kh  -  kh  cos  kh)  d/X  <  m/N,  (19) 

K 

T'(m)/n\ 

p)'")  ~  2(1  -  cos  kh)  d/X  <  m/N.  (20) 

n. 

so  that  and  Pj”‘^  arc  slowly  varying  functions  of  kh  that  also  possess  the 

Properties  (2)  and  (3)  above. 

C.  Phase  Sequence  Resonances 

Throughout  this  section  d/X  is  the  variable  and  a/X.  h/X  and  .V  are  fi.xed,  just 
as  in  the  figures  in  [G].  The  physical  picture  is  similar  when  the  frequency  is  varied 
and  the  geometrical  parameters  a,  h.  and  d  are  fi.xe<l. 
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A  typical  plot  of  +  as  afunction  of  d/X  is  given  in  Fig.  1.  It  is 

seen  that  is  a  quantity  of  order  1  that  has  two  zeros  in  the  range  0  <  d/ X  <  0.5. 
The  array  is  at  its  phase- sequence  resonance  when  is  exactly  zero.  In 

Fig.  1,  the  smaller  root  is  not  in  the  region  of  validity  of  the  two-term  theory  since 
d  <  h.  Since  the  position  of  the  resonance  is  determined  only  by  the  real  part  of  the 
kernel,  a  particular  m/N  phase-sequence  resonance  will  occur  at  roughly  the  same 
value  of  d/X  for  all  large  N.  Hence,  it  is  meaningful  to  examine  resonant  currents, 
etc.,  as  N  becomes  larger,  while  keeping  m/N  and  d/X  fixed.  Denote  by  the 

position  of  the  larger  root,  so  that  =  0  when  d/X  =  It  is  seen  from  (3) 

and  (15)  that,  at  the  phase-sequence  resonance, 


)(">) 


d('«) 
R _ 


>("*) 


D\ 


^  res 

= 


2ir 


pirn) 

_ ^  R 

Co'^dRCoskh  /)('«)’ 

27r 


p(m) 

sin  kh  +  (1  —  cos  kh)  — - — 


D 


(m) 


(21) 

(22) 

(23) 


where  Pr^\  and  aie  evaluated  at  d/X  = 

The  quantity  P^^^  is  of  order  1.  Because  of  Properties  (2)  and  (3)  and  equations 
(19),  (20)  and  (16),  if  <  rn/N,  it  is  seen  that 


Property  (4):  At  the  m^*'  phase-sequence  resonance,  the  phase- 
sequence  conductance  is  exponentially  large 

in  N .  The  self-  and  mutual  conductances  Gj  ^ 
around  the  array  are  also  exponentially  large 
and  they  vary  around  the  array  according  to 


Gj  /  a  cos 


27r(/  —  l)/n 
N 


;  /  =  1,  2,  ...,  iV. 


(24) 


Property  (5): 


GIcs*  and  the  Gj  /s  will  be  mvich  larger  when  the 
resonant  spacing  d/X  —  at  a  smaller  value. 
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On  the  other  hand,  and  the  /s  are  not  large  when  the  array  is  exactly  at 

resonance.  The  magnitudes  jFj  ,|  of  the  self-  and  mutual  admittances  are  determined 
entirely  by  the  conductances,  so  that  Kj  /  oc  |cos[27r(/  —  l)m/A'^]|.  The  phases 
Arg(yj  j)  can  be  either  0°  or  1S0“;  in  the  special  case  when  cos[27r(/  —  l)m/iV]  =  0, 
the  current  on  element  I  is  very  small  compared  to  that  on  all  other  elements.  These 
observations  are  generally  consistent  with  the  figures  in  [6].  The  differences  are  due 
to  the  fact  that  in  [6]  the  resonances  are  not  as  narrow. 

The  parameter  depends  on  a/\  only  through  and 

on  d/X  only  through  D^j^.  is  zero  when  the  elements  are  self- resonant.  By 

plotting  for  various  values  of  /i/A,  it  is  seen  that  is  a  decreasing  function 
of  a/A,  at  least  when  0.001  <  a/A  <  0.07  and  a/A  <C  h/X  <  0.22.  (It  can  be 
shown,  in  fact,  following  a  procedure  similar  to  that  in  Appendix  II  of  [3],  that  the 
variation  with  a  is  linear  in  the  quantity  Q  =  21n(2h/a},  but  the  approximation  is 
poor  when  h/a  becomes  small.)  Hence,  making  the  dipoles  electrically  thicker  will 
result  in  decreasing  the  anrplilude  of  a  curve  like  that  in  Fig.  1,  thereby  shifting 
the  resonant  spacing  to  a  smaller  value  of  d/X.  The  resonant  currents  will 

therefore  become  much  larger.  When  the  elements  are  electrically  very  thin,  the 
array  can  have  no  narrow  resonances  at  all,  because  the  resonant  root  occurs  at  a 
value  d/X  >  m/N .  This  is  the  case  in  Fig.  1  of  [6]. 

The  effect  of  changing  the  length  h/X  is  much  more  involved,  since  both 
and  depend  on  h/X  in  a  complicated  way.  However,  extensive  numerical 

calculations  show  that  the  position  of  the  root  decreases  when  h/X  increases, 

at  least  when  a/A  and  h/X  are  in  the  above-mentioned  ranges.  Table  1  shows 
the  resonant  spacings  the  values  of  K\^\0)  evaluated  at  d/X  =  and 

the  self-conductance  Gj  j,  for  90-elcment  arrays  at  their  ?n  =  N/2  phase- sequence 
resonance  as  a/A  and  h/X  vary. 


Table  1 


Resonant  Spacings  d/X  = 

Values  of  the  Imaginary  Part  of  the  Kernel  at  z  =  0,  djX  = 

and  the  Driving- Point  Conductance  Gj  j  at  Resonance. 
Number  of  Elements  N  —  90;  Phase  Sequence  m  =  N/'2  =  45. 
Roots  are  Sought  in  the  Intei'val  h/X  <  d/X  <  X  —  jTt/J\\ 


h/X 

a/X  =  0.01 

a/X  =  0.03 

a/X  =  0.05 

no  root 

no  root 

=  0.479 

0.14 

=  -0.25 

A 

Gj  1  =  16.6  mA/V 

no  root 

^m/N  =  0-4^^0 

^m/N  =  0-439 

0.16 

= -0.1S 
k 

=  -4.8  X  10-3 

Gi  1  =  11.8  mA/V 

Gi,i  =  132  mA/V 

K/n  = 

K/n  =  OWl 

K,/n  =  0-370 

0.18 

k 

= -2.1  X 10-’ 

’“''y  =  -3.9  X  10-3 
k 

Gi  i  =  5.2  mA/V 

Gj  1  =  196  mA/V 

Gi,j  =  9.6  X  10^  mA/V 

K/n  =  0«7 

^m/N  —  0-336 

■»,„/«  =  0-273 

0.20 

■  ^  ,  ’  =  -4.1  X  10"^ 

k 

=  -7.0  X  10-.» 

K 

=  -2.3  X  10-« 
k 

Gi^i  =  90  mA/V 

Gi  1  =4.2  X  10»  mA/V 

Gj  1  =  1.2  X  10^5  mA/V 
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In  Table  1,  the  imaginary  part  of  the  kernel  was  evaluated  from  (12)  using 
quadruple  precision  and  the  self-conductance  was  calculated  using  the  exact  formu¬ 
las  (1)-(16),  with  set  equal  to  zero.  The  extremely  large  currents  when  the 

elements  are  long  and  thick  require  perfectly  conducting  elements;  they  could  be 
realized  in  practice  only  with  superconducting  elements.  The  conclusion  is  that 


Property  (6):  If  a  specific  phase-sequence  resonance  is  desired, 
making  the  dipoles  electrically  longer  or  thicker 
will  require  an  electrically  smaller  circle  and  will 
result  in  much  higher  resonant  currents  around  the 
array,  at  least  when  a/A  and  h/X  are  in  the  ranges 
0.001  <  a/A  <  0.07  and  a/A  <  h/X  <  0.22. 


Of  course,  the  resonant  currents  have  a  complicated  explicit  dependence  on  h/X 
and  a/A  as  well  as  this  implicit  dependence  through  d/A  =  but  the  latter  is 

dominzuit. 

D.  Behavior  near  a  Phase-Sequence  Resonance 

Consider  again  that  a/A,  h/X  and  N  are  fixed  and  that  d/X  is  varied,  but  very 
close  to  a  resonant  spacing  array  is  very  close  to  its  m‘**  resonance. 

The  function  is  usually  a  quantity  of  order  1,  but  near  resonance  it  is  of  the 

order  of  magnitude  of  the  exponentially  small  quantity  it  is  the  controlling 

quantity  in  (3).  It  is  a  good  approximation  to  assume  that  P\^^  and 

are  constant  and  that  '  varies  linearly  with  d/A;  ’  oc  d/A  —  Using 

^  Pj"'\  it  is  seen  from  (3)  that  Re{T^'"^}  has  extrema  when 
with  the  corresponding  values: 

p(”>)  1 

RefT*™*}  =  ±-^  =  (25) 

2D  i  “ 


From  the  relations  (15)  between  T^'^\  and  and  (16)  between  j  and 
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B^'^\  it  is  seen  that 

Property  (7):  and  the  B^  ,’s  are  very  rapidly  varying  near  a 

narrow  resonance.  When  the  spacing  d/X  is  such 
that  has  decreased  to  half  its  maximum  value, 

is  equal  to  =  ±|G^^s^ 

Hence,  B^'”^  and  the  Bj  /s  have  a  zero  very  close 
to  resonance.  The  Bj  ^’s  vary  around  the  array  as 
cos[27r(/  —  l)m/A^]. 


It  is  therefore  possible  to  design  an  array  near  resonance  with  a  purely  resistive 
driving-point  impedance,  but  this  property  will  be  extremely  sensitive  to  slight 
changes  in  the  parameters.  It  is  interesting  to  note  that  the  circular  array  is  similar 
in  behavior  to  a  single  center-driven  dipole  [13].  In  both  cases,  when  the  driving- 
point  conductance  is  maximum  as  a  function  of  /,  the  susceptance  is  nearly  zero. 
When  the  element  radius  becomes  larger,  the  Q  rises. 

The  Q  of  the  resonant  array  was  defined  in  [6].  It  may  be  estimated  from  the 
curve  of  Gj  ^  as  a  function  of  d/X  as 


Qr 


^m/N 

id/X),-{d/X),- 


(26) 


(The  actual  definition  involves  the  frequency.)  In  (26),  the  {d/X)2  and  (d/A)j  are 
the  spacings  at  which  the  power  is  reduced  to  one-half  the  maximum  at  constant 
voltage,  i.e.,  when  Approximating  as  a  linear  function,  viz., 

~  a{d/X  —  using  ^  d\’^\  it  is  seen  that 


Qr  ^ 


(m) 


(27) 


(£)(m))2- 

It  should  be  noted  that  it  is  very  difficult  numerically  to  calculate  quantities  near 
resonance.  is  the  difference  between  two  nearly  equal,  complicated  integrals 

of  order  1  and  must  be  evaluated  with  very  high  precision.  Formula  (20)  provides 
a  simple  way  to  estimate  the  Q.  A  formula  similar  to  (20)  for  the  bandwidth  may 
also  be  obtained. 
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E.  Radiation  Field  at  or  near  Resonance 

With  the  two-term  theory  currents,  and  with  the  center  of  the  spherical  co¬ 
ordinates  (r,6,(f))  placed  at  the  center  of  the  array,  the  radiation  field  is  given  by 
[11]: 


Eir,e,<f>)  =  9Es 


'I'j^cosA:/i  r 


X  i  J'^gjtpsin  ecos((/.-<^,)  I 

^  /=1  ^ 


where  (f>i  =  27r(/  —  1)/N  is  the  location  of  element  /,  and 


G(e)  = 


cos(kh  cos  6)  —  cos  kh 
sin  9 

sin  kk  cos(kh  cos  9)  cos  9  —  cos  kh  sm{kh  cos  9) 
sin  9  cos  9 


In  (28),  the  first  term  represents  radiation  from  the  sine  current  of  the  driven 
element;  the  second  term  is  radiation  from  the  shifted-cosine  currents  of  all  elements. 
At  or  near  a  narrow  resonance,  we  have  Tf  ~  Tj  cos[27r(/  —  \)m/N].  It  will  be  seen 
that  the  first  term  in  (28)  may  be  neglected,  so  that  it  is  convenient  to  define  the 
array  factor  for  the  m‘''  phase-sequence  resonance  as 


1  ^ 
i^i 


~"(^  jkpsin  0cos(<t>-<i>,) 

cos  t 


(31a) 


(31b) 


so  that  the  radiation  field  is  approximately 


=  -,T,  II  —G(e)T,A'’"\0,4:). 

cos  Ah  r 


(32) 
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Thus,  the  array  factor  is  the  radiation  field  due  to  a  circular  array  of  isotropic  radi¬ 
ators  with  the  phase-sequence  resonance  currents  around  the  array,  element  1 
having  unit  current.  The  sum  in  (31b)  may  be  written  exactly  as  follows;  see  [14] 
for  a  detailed  derivation: 


=  Nj^JJ^N{d/X)s\n6]  cos(m^6) 

OO 

+  iV  J^p_,„[iV(d/A)sin^]  cos[(A>  -  m)cf>] 

P=1 

OO 

+  N  53 >'''■+"•  J;vp+mlA'(<</A)sin«l  cosKA'p  +  m)^].  (33) 

p=l 

Because  of  the  condition  d/X  <  m/N  <  the  arguments  of  the  Bessel  functions 
are  always  smaller  than  the  orders.  When  N  is  large,  only  two  terms  in  (33)  are 
significant,  viz., 

~  Nj^J^[N(d/X)sin6]  cos{m<f>) 

+  {same  with  rn  —*  N  —  m}.  (34) 

As  with  the  imaginary  part  of  the  kernel,  the  first  term  is  adequate  when  m  <C  N /2 
and  the  second  term  is  equal  to  the  first  one  when  m  =  N/2.  Assuming  for  simplicity 
that  m  N/2  and  using  the  asymptotic  formula  for  the  Bessel  functions,  one 
obtains 


(2irm/.V)'/2  (l-(i,„sin«)21'/< 

X  exp[— A^(m/A^)<7(x„,  sin^)]  cos(m<^). 


where  and  g(x)  are  defined  in  (18).  Hence,  the  array  factor  is  an  exponentially 
small  quantity  and,  in  fact,  it  shares  Properties  (2)  and  (3)  of  the  imaginary  part  of 
the  kernel.  This  is  due  to  cancellation.  Also,  it  has  a  zero  of  order  rn  at  6  =  0.  From 
(35),  the  field  formula  (32),  and  the  expressions  (19)-(21)  for  Tj  ~  {2/N)Tr^\  we 
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see  that: 

Property  (8):  The  magnitude  of  the  radiation  field  at  any  fixed 
point  in  space  is  exponentially  large  in  N. 


That  is,  the  resonant  currents  are  large  enough  to  overcome  the  cancellation  effects. 
Radiation  from  the  sine  current  of  element  1  is  negligible. 

Property  (9):  The  horizontal  field  pattern  {9  =  7r/2)  consists  of 
2m  spikes. 

Property  (10):  The  vertical  field  pattern  is  very  narrow,  with 
a  maximum  at  ^  =  7r/2. 


An  exponential  largeness  of  the  field  should  be  expected  since,  for  lossless  ele¬ 
ments,  integration  of  over  a  large  sphere  should  give  the  total  radiated  power 
which  is  exponentially  large.  The  narrowness  of  the  vertical  beam  can 
be  estimated  by  neglecting  variations  of  the  field  in  (32)  due  to  the  slowly  varying 
G{9)  and  defining  vertical  directivity  as  the  maximum  of  the  array  factor  divided 
by  its  mean  value,  viz.. 

Subject  to  the  approximation  (34),  the  integral  can  be  carried  out  analytically. 
Thus, 


„ _ •2J„\N(dlX)\ _  2J„[JV(rf/A)l 

where  the  leist  approximation  is  valid  for  large  m.  With  the  asymptotic  expression 
for  the  Bessel  functions. 


Dy  ~  (2iV/7r)'/2[(m/iV)2  -  (d/A)^]*/'*.  (38) 

Hence, 

Property  (11):  For  a  specific  ph^lse-sequence  resonance  m/N, 
making  N  larger  will  result  in  a  narrower 
vertical  field  pattern,  and  in  more  spikes 
in  the  horizontal  plane. 
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Property  (12):  For  fixed  N  and  for  a  specific  phase- sequence 
resonance  m,  making  the  dipoles  thicker  or 
longer  will  result  in  a  smaller  resonant  spacing 
^m/N’  ^  much  narrower  resonance,  and  a 
slightly  more  directive  vertical  field  pattern. 

The  directivity  may  therefore  be  made  arbitrarily  large  by  making  N  large  (although 
the  increase  is  slow,  as  the  square  root  of  N).  The  field  strength  at  any  point  in  space 
increases  exponentially.  The  input  impedance  may  be  a  pure  resistance.  How('ver. 
the  physical  dimensions  of  the  array  increase  (linearly  with  N)  and  the  bandwidth 
decreases  exponentially. 

The  array  factor’s  smallness  has  an  interesting  consequence.  For  resonant  non¬ 
circular  arrays,  an  array  factor  A{d,4>)  may  be  defined  exactly  as  in  (31a).  This 
array  factor  will  depend  on  the  array’s  geometry  and  the  relative  current  distribu¬ 
tion  around  the  array.  It  will  be  a  sum  of  N  terms  of  order  1,  each  term  depending 
on  the  location  of  element  I  and  its  relative  current  (admittance).  If  the  noncircular 
array  is  thought  of  as  a  perturbation  of  some  corresponding  circular  array,  then 
it  is  logical  to  assume  that  the  current  distribution  around  the  array  will  not  be 
significantly  affected.  Hence,  each  term  in  the  sum  for  A(9,  4>)  will  be  close  to  each 
term  in  the  sum  for  the  circular  array.  However,  any  very  small  quantity  that  can 
be  written  as  the  sum  of  terms  of  order  1  is  extremely  sensitive  to  perturbations 
of  these  terms.  Therefore,  the  array  factor  (field  pattern)  for  the  nonrircular  array 
will  not  be  close  to  that  of  the  circular  array.  The  conclusion  is  that  any  pertur¬ 
bation  of  the  resonant  circular  array’s  shape  resulting  in  only  a  slight  change  in 
the  current  distribution  around  the  array  will  significantly  change  the  field  pattern. 
This  means  that  a  wide  variety  of  field  patterns  may  be  obtained  by  resonant  non¬ 
circular  arrays,  perhaps  even  a  superdirective  pattern.  However,  one  must  be  very 
careful  in  dealing  with  noncircular  arrays.  The  N  linear  algebraic  equations  for  the 
currents,  to  which  the  N  coupled  integral  equations  reduce  by  the  two-term  theory 
approximations,  must  be  solved  with  high  accuracy.  The  currents  must  be  known 
precisely  in  order  to  obtain  any  reasonable  estimate  of  the  far-field  pattern. 
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